Quantum coherent control of highly multipartite continuous-variable entangled states 

by tailoring parametric interactions 
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The generation of continuous- variable multipartite entangled states is important for several proto- 
cols of quantum information processing and communication, such as one-way quantum computation 
or controlled dense coding. In this article we theoretically show that multimode optical parametric 
oscillators can produce a great variety of such states by an appropriate control of the parametric 
interaction, what we accomplish by tailoring either the spatio-temporal shape of the pump, or the ge- 
ometry of the nonlinear medium. Specific examples involving currently available optical parametric 
oscillators are given, hence showing that our ideas are within reach of present technology. 

PACS numbers: 42.50.Dv-Quantum state engineering and measurement, 42.65.Yj-OpticaI parametric oscil- 
lators and amplifiers, 42.65.Re-Ultrafast processes, optical pulse generation and pulse compression. 
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I. INTRODUCTION 

Squeezed states of light were introduced several 
decades ago as states which could help beating the limits 
set by quantum mechanics on the precision of measure- 
ments performed with coherent light [11, [1]. These are 
states in which one of the quadratures of light (equiva- 
lent to the position and momentum of a mechanical os- 
cillator) has a quantum uncertainty below the vacuum or 
"shot noise" level, at the expense of increasing the quan- 
tum fluctuations in the orthogonal quadrature. These 
states found their way into the new century thanks to 
their widely proved applications in ultra-precise metrol- 
ogy (such as gravitational wave detection [3, 0] or beam 
displacements [1, Q ) , as well as in the field of quantum 
information with continuous variables 0, where the 
highest-quality entangled states (the basic ingredient of 
many quantum information protocols) known to date are 
currently obtained by mixing squeezed beams with linear 
optics 

Squeezed light can be obtained via the parametric 
down-conversion process that takes place inside a sec- 
ond order nonlinear medium pumped by a laser beam 
In order to increase the nonlinear interaction, it is 
customary to insert the nonlinear crystal in an optical 
cavity — dealing then with a so-called "optical paramet- 
ric oscillator" (OPO) — , and large levels of squeezing are 
obtained in the down-converted field when the OPO is 
operated close to threshold, 93% of noise reduction be- 
ing the current benchmark [ll| (see also [l2l. ITsj). 

On the other hand, quantum information has reached a 
stage where real-world applications stimulate an intense 
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research for the implementation of reliable and practical 
quantum protocols for quantum communication and in- 
formation processing. Several of the promised benefits 
require though a quantum "substrate" that is created by 
distributing quantum correlations (entanglement) among 
a number of degrees of freedom (modes) increasing with 
the complexity of the task to achieve. The protocols of 
quantum telecloning [13, HBl and controlled dense coding 
[igI [iTj constitute paradigmatic examples of this scenario 
for small number of modes, while one-way quantum com- 
putation |18 21], in which the computation is achieved 
by applying local measurements to a set of modes ini- 
tially in a cluster state, is a most promising example in 
the large number of modes regime. 

However, the generation of such multipartite entangled 
states by means of optical devices requires experimental 
configurations whose cornplexity increases with the num- 
ber of modes involved [l3, [13 [22l - [25j . In contrast, a 
practical source should be compact, scalable, and permit 
to master the quantum properties of the generated states 
even when the number of modes is very large. 

A continuous-variable cluster state source with these 
properties was proposed in [2^ — and refined in subse- 
quent papers [13, [11] (see also [l^]) — , which consisted 
in an OPO driven by a multifrequency pump field in 
the presence of concurrent nonlinearities. Recently, these 
ideas have been brought to the laboratory, and the record 
of 15 quadripartite cluster states have been generated 
aloiig 60 consecutive longitudinal modes of the OPO cav- 
ity [23 |. This interesting approach is however limited by 
the number of frequencies that can be coupled within 
the phase-matching bandwidth, the need of frequency- 
sensitive measurements, and the fact that a nonlinear 
crystal with different phase-matching properties has to 
be engineered for each particular cluster state that one 
desires to create. An even more simple and promising 
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approach based on the synchronously operation of only 
one single mode vacuum squeezer and a quantum non- 
demolition gate was proposed very recently in [soj . 

In this paper we propose an alternative approach based 
on the use of a naturally multimode optical parametric 



oscillator, either in the spatial |31| - |42| or in the temporal 
domain [i^, . Below threshold, they produce multi- 
mode squeezed states, which has been observed experi- 
mentally [l^jllll- When operated above threshold, such 
multimode devices have been recently proved to be ca- 
pable of producing non-critically squeezed states of light 
via the phenomena of spontaneous symmetry breaking 
[H m [si EH, 13 and pump clamping [H El . 

The advantage of the scheme introduced in the current 
article for the generation of cluster states over previous 
proposals [20l . [26l - l28l | is that one does not need to engineer 
the couplings between each mode through crystal design 
(which is specific of each coupling scheme) : choosing in- 
stead the temporal or spatial shape of the pump (which 
is a flexible technique) is enough to directly control the 
characteristics of the generated quantum state, even for 
a very large number of modes. We will show that it is 
possible to master in this way the number of entangled 
modes, as well as the distribution of their quantum corre- 
lations, thus permitting the generation of arbitrary mul- 
timode Gaussian quantum states of any dimension. The 
technique that we propose is therefore reminiscent of the 
widely used "coherent control" of evolution of atoms and 
molecules that is obtained by appropriately shaping the 
pulses with which they interact [43, |48| . 

The article is organized as follows. In Section |ll] 
we show how both the light generated via single-pass 
spontaneous parametric down conversion in the low-gain 
regime, as well as that leaving a multimode optical para- 
metric oscillator below threshold, can be described as 
a combination of independent, squeezed modes (termed 
"supcrmodes" later on). After considering a particular 
(but physically relevant) example of the parametric inter- 
action to show how these supermodes arise (Section llll| . 
we demonstrate that the squeezing and spatio-temporal 
shape of such modes can be controlled either by tailoring 
the shape of the pumping field, or the geometry of the 
nonlinear medium; we consider separately the cases of an 
OPO with either many longitudinal modes fSection llV A[) 
or just a few transverse modes (Section IIVB|) available 
for parametric down-conversion. We then explain how 
our results can be used for the generation of arbitrary 
cluster states of copropagating modes in Section |Vl and 
give our conclusions in Section IVll 



II. GENERALIZED SUPERMODE 
DESCRIPTION OF WEAK PARAMETRIC 
INTERACTIONS 



nel function or matrix (depending on the continuum or 
discrete nature of the modes involved) which describes 
the coupling among the different relevant signal-idler 
modes. As shown in previous studies, the diagonaliza- 
tion of such kernel is instrumental for determining the 
objects with well-defined quantum properties, which turn 
out to be linear combinations of the signal-idler modes. 
Such linear combinations were termed "supcrmodes" in 
[43| and the next subsections are devoted to their intro- 
duction in both single-pass and cavity devices. We will 
consider collinear type I degenerate phase matching for 
definiteness but our treatment can be generalized easily 
to other types of phase matching. 



A. Single-pass parametric interactions 

In the single-pass configuration the pump beam ampli- 
fies paranictrically the quantum noise impinging the crys- 
tal around the subharmonic frequencies and the down- 
converted field is measured at the exit of that crystal. 
In the undepleted pump approximation, holding when 
parametric gain is low, the "output" boson operators are 
generically related to the "input" ones via the well known 
Bogoliubov transformation 



aout(C) = / dC' \C (e, C) ai„ {O + S it a al (^') 



(1) 

where ^ is a general continuous index standing for fre- 
quency or/and transverse position^, and C and S are 
Green functions that solve the propagation equations 
(see [13 for details). In the very weak conversion limit 
a perturbative treatment can be applied to first order 
in the crystal length, which physically means consider- 
ing only generation of single photon pairs. In such case 
C (e, C) = S{^- e') and S (^, = gK ($, C), where g is 
a coupling constant proportional to the length Ic and non- 
linear susceptibility x^^^ of the crystal and to the square 
root of the total pump irradiance P, and A'(^,^') is a 
kernel. Hence in this regime ^ becomes 



aout(0 « ai„(0 + 9 / dCKi^,a&l iO 



(2) 



As we will see immediately, the form of K is universal, 
and is given by the product of the pump amplitude and 
some function describing the overlap between modes over 
the crystal. However, the special form of K depends on 
the considered problem and in the following we treat suc- 
cessively the temporal and spatial aspects of single-pass 
parametric interaction in the single photon pair approx- 
imation ([2]). 



The dynamics of parametric interactions in the low- 
gain regime (for single-pass devices) or below its oscilla- 
tion threshold (for cavity devices) is controlled by a ker- 



^ Transverse means orthogonal with respect to a propagation di- 
rection, in which case the paraxial approximation is assumed. 
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Let us begin by the spectral/temporal aspects in which 
case ^ ~ oj corresponds to the frequency of the monochro- 
matic modes (we call this the temporal case) . Assuming 
that all the parametrically coupled modes have the same 
transverse dependence, and neglecting diffraction inside 
the crystal, the coupling kernel K(uj,uj') turns out to be 



K{uj, uj') = ap (w + uj') D {uj, uj') , 



(3) 



Q!p (w) being proportional to the spectral pump ampli- 
tude at frequency w, and D the usual phase- matching 
function^ 



D{uj,oj') 



/e/2 



a;, a; ]z 



sine [$ {uj,uj')] , 



-lc/2 



(4) 

where Afc (w, uj') = k {uj + uj') — k (uj) — k {uj'), k (uj) is 
the optical wavenumber at frequency uj inside the crystal, 
sine (x) ~ sin (x) /x is the sinus cardinal function, and 



$(cj,cj') = Ak{uj,uj')lc/2, 



(5) 



is a phase mismatch. 

Let us now turn to the spatial case within the paraxial 
approximation, in which case ^ = r is the transverse 
spatial variable (d^ = d^r). We define this "spatial case" 
in the sense that pump is assumed monochromatic and 
focus is put on the signal/idler (multimode) field at just 
the subharmonic frequency, in which case the kernel is 
given by [40| : 



if(r,r') = ap(^^)A(r-r- 



(6) 



ttp (r) being the normalized pump amplitude at trans- 
verse point r, and A(r) the diffraction function 



A(r) = 



ifcs , 


dz { 
— exp 


ifcs 1 


1 




r 


■)1 




^coh 





(7) 



ks being the phase-matched signal wavenumber inside 
the crystal, Zcoh = \/2Zc/fcs the "coherence length" , and 
Si {z) = sine (u) du the sine integral function. 

In the general case, the spatial and spectral aspects 
of the parametric interaction are simultaneously present, 
giving rise to new interesting features (50| which we will 
not consider here. 



B. Intracavity parametric interactions 

When an optical cavity is used to enhance the efficiency 
of the nonlinear process, the previous nonlinear couplings 
represented by kernel K are projected onto the cavity 
modes, and the multi-dimensional spectral and spatial 
properties of the correlated photons can be lost due to 
such filtering. However for a given length and geometry 
of the cavity, which we assume to have a cylindrical sym- 
metry around the optical axis, a great number of modes 
can be simultaneously sustained: in the frequency do- 
main they are the series of longitudinal modes separated 
by the free spectral range of the cavity; in the transverse 
spatial domain they are the set of Laguerre-Gauss modes 
{TEM,,}f,f. 

In the temporal case we assume the pump consisting 
of an unlimited series of pulses at a given repetition rate, 
equal to the free spectral range of the cavity: this is 
what is called a Synchronously Pumped OPO (SPOPO) 
[4^ 13]. The pump spectrum is thus a frequency comb, 
consisting of a large number of frequency components, 
each of which gives rise to signal-idler photons belonging 
to different longitudinal modes through the parametric 
down conversion process, and selected by energy and lin- 
ear momentum conservation. 

In the spatial case, on the contrary, the pump is as- 
sumed monochromatic with a given spatial profile. In 
this case, owing to the linearity of the interaction in the 
below threshold regime, one can focus on, say, the sub- 
harmonic signal photons. However these photons belong, 
in general, to different transverse modes TEMpj (degen- 
erate in frequency) , characterized by a constant value of 
the sum / = 2p + | / 1 , known as "family index" [s^l • As 
well, even several families of transverse modes can be- 
come relevant in the case of a degenerate cavity such as 
the the confocal [3l|, [13, [5lj or the self-imaging cavity 
[40I [5^ , in which many different families (with different 
family indices /) resonate at the same frequency. We will 
use a single generic index to label this discrete series of 
modes. 

It is instructive in this intracavity interaction prob- 
lem to write down the interaction Hamiltonian in the 
undepleted pump approximation, which can be written 
in general as 



Hi^i^ I dCdC'/v it a iO in + H.c, 



(8) 



where the continuous boson operators a(^) satisfy the 
standard commutation relations 

[m,m] = [«no,«nc')] = o, (9a) 



^ In [49I a factor exp {i [k (lj) — k (oj')] Zc/2} is included in the ex- 
pression Here it is absent because we are implicitly working 
in the interaction picture with respect to time and space, analo- 
gous to Eq. (17) in see [11,13 for details. 



and the kernel K (^, ^') plays a role analogous to that in 
single-pass devices. 

In the temporal/spectral case (SPOPO) the Hamilto- 
nian of the nonlinear interaction can be written in terms 
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of discrete creation and annihilation operators of the cav- 
ity longitudinal modes: 



(10) 



where Kij are the specialization of the continuous- 
variable kernel K ([3]-[5]) to the relevant longitudinal 
modes around the subharmonic, and cii is the boson op- 
erator for the signal- idler mode i verifying [a^, a^] = 6i^j. 
The dynamics of the intracavity modes is then given by 
the following quantum Langevin equations 



doi 
"d^ 



'lai 



7a 



(11) 



where 7 is the cavity loss rate (assumed identical for all 
signal/idler modes - which, by the way, cannot be distin- 
guished from one another in this coUinear tjpe I degen- 
erate phase matching case), and cr^ = (5/7) is a dimen- 
sionless pumping parameter proportional to the actual 
pump power (see, e.g. [i^l for more details). 

In the spatial case (as defined in the previous Section) , 
under circumstances in which only a few signal modes are 
relevant [s^, Hamiltonian ([TU]) holds by identifying the 
indices («,j) with the available TEMpj modes at the sig- 
nal frequency, whose dynamics are ruled by the quantum 
Langevin equations (|lip . In degenerate cavities, such as 
confocal or self-imaging cavities, a discrete representa- 
tion is still possible, but a continuous one ([8]) is more 
helpful, in which case the quantum Langevin equations 
become 14(1 



da{T) 
dt 



--fa (r) + V27ain (r)+7cr J d^r'A' (r, r') (r') • 
1^ 

For self-imaging cavities l40j the kernel is given by (|6]) 
and ([7]), while for the confocal case it has a slightly dif- 
ferent expression |3l| owed to the fact that even and odd 
transverse families resonate at different frequencies. 



C. Supermodes 

We have seen in the previous Section that the output 
from parametric devices in the low-gain (or below thresh- 
old) regime is governed by a kernel K, which is equal, 
under equivalent conditions, in the single-pass and in- 
tracavity cases. As both the input-output relations in 
single-pass configurations and quantum Langevin equa- 
tions in cavity devices are linear in the considered regime, 
diagonalization of the kernel allows a considerable sim- 
plification in the description of the problem as well as 
a clear physical picture of the entities in which clean 
quantum properties are concentrated. This approach 
pol I43I Hj, [4^ is the single-partite and continuous- 
variable version [53| of the Schmidt decomposition for the 
bi-partite bi-photon wave function [s^ , which in general 



is formally expressed by the Bloch-Messiah decomposi- 
tion m. 

Let us introduce the eigenmodes of the continuous ker- 
nel i^(^,^'), which are the solutions Sn{0 of Ih'^ Fred- 
holm integral equation: 

F [s„ (e)] = / dc if (e, r)s„(e') = a„ ,s„(e). (la) 



Starting from the eigenmodes s„, we can define the 
associated supermode annihilation operator Sn as 



5„ = J d^sniom, 



(14) 



where we assume Sn(C) to be normalized so that 
[Sn, S^] = 6n.m- Thcsc Operators annihilate a photon 
in a mode which is a combination of either spatial or fre- 
quency modes, and we call it "supermode" [i^ because 
it is the combination of many cavity modes and has a 
non trivial spatial or spectral/temporal shape in general. 
Using (dH) and ([131) into Eq. © we get 



So 



Su 



- gAnS' 



t 

in.n ■ 



(15) 



When the crystal is inserted in an optical cavity — OPO 
case — , Ref. [4j] shows that the supermodes can be de- 
fined as a discrete combination of the cavity eigenmodes 
instead of the continuous combination written above, cor- 
responding to the diagonalization of matrix K in Eq. 
([TT|). These two possible expressions are identical in the 
limit where the cavity modes are tight enough so that the 
sum can be assimilated to an integral. In the following 
we assume that this condition is fulfilled. In this case, the 
evolution of these supermode operators is governed by a 
set of uncoupled quantum Langevin equations obtained 
by introducing the discrete versions of ([H)) and in 
Eq. (HH): 



dSn A 
-df = 



2-rSn 



■7 



TAnSl 



(16) 



This equation shows that according to this definition, 
any multimode OPO can be in fact seen as a set of inde- 
pendent single-mode OPOs, and therefore can produce a 
set of copropagating squeezed supermodes. Hence each 
supermode is independently squeezed, and the squeezed 
quadrature can be shown to have a noise level at zero 
frequency equal to (|Ai| - |A„|)2/(|Ai| 4- |A„|)2 (1 setting 
the shot-noise level) when the system is operated close to 
threshold, where Ai is the eigenvalue of largest absolute 
value. 



III. A GENERAL METHOD FOR TAILORING 
THE SUPERMODE SPECTRUM 

The solution of Eq. is simple when the kernel 

factorizes, 



(17) 
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in which case the unique cigcnmode is, within a mul- 
tiphcative factor, si = /, with corresponding eigen- 
value Ai = /d^' /($') ■ This resuh can be extended to 
the case where there is a set of N orthogonal functions 
{fn}n=i such that 

N 
n=l 

in which case Eq. (|13l) has N solutions .s„ = /„ (save 
multiplicative constant), A„ = / d^'/„ (^') . When K 
has the same type of decomposition as (jlSp . but with 
non orthogonal functions, just linearly independent, one 
can show that (|T3)) has still N solutions [5g, which are 
now different from the functions /„. 

The analytical description of the kernel K in terms of a 
basis of linearly independent functions (supermodes) like 
in Eq. (|18p has a very simple solution and physical inter- 
pretation in the case where the kernel can be factorizcd 
as 

K{x,x') = K+{x + x')K^{x- x'). (19) 

In the spatial case this is the form encountered [13], see 
dH]): K+ relates to the pump, while K- has to do with 
the crystal. In the temporal case, while this is not the 
most general kernel one can find in actual applications, 
it has been shown in [13, that it is a sensible approx- 
imation to many real cases. The function in (fT9|) 
can be manipulated by tailoring the pump shape, while 
in order to tunc the function A'_ one needs to play with 
the geometry of the nonlinear crystal. 

We now exhibit a series of "simple" kernels which allow 
a high degree of control over their spectra and, as wc 
will show in the next Section, find applications in actual 
systems. 

Let us start by the simplest case, the symmetric Gaus- 
sian kernel, defined by K+{x) ~ A'_(x) = e~2°' ^ , read- 
ing 

Kix,x') =c-'''''\~''^^'\ (20) 

as in ()17p . whose only supermode is the Gaussian func- 
tion si (x) — e~'^ ^ , with corresponding eigenvalue Ai = 

Let us now consider the previous symmetric kernel 
(|20p . now multiplied by cosine functions of x + x' and 
X — x\ i.e.: 

JV± 

K± ix) = e-i'^'-' cos {f3tx) , (21) 

n=0 

which is a quasiperiodic modulation of the Gaussian ker- 
nel. In ((2T|) . {bn,Pn}n=o ^'^^ coustauts and = by 
definition. This leads to a kernel K ([T9|) which can be 
written like in Eq. (|18p — with non-orthogonal functions 
/„ — using well known trigonometric formulae. Wc have 



therefore shown that a kernel having the form ([19]) with 
(f2T|) win have exactly (2A^+ -|- 1)(27V_ -I- 1) eigcnmodcs, or 
AN^N- modes if the constant terms = 0. As shown 
in Appendix A a simple way to diagonalize the kernel 
defined by (|2T|) consists in using auxiliary functions 

4V,*4 (^) = c-'^'-'ti (/?+ x) h {p-,x) , (22) 

where ti=i,2 stands for any of the trigonometric func- 
tions cos or sin, in terms of which actual eigenvectors 
and eigenvalues can be found by linear combinations of 
them. The expressions for eigenvalues and eigenvectors 
are awfully cumbersome but one can show that whenever 

{(3:^y,{i3-y»8a^ (23) 

the functions (j22[) are approximate eigenvectors indeed. 
For each couple (ni,n2) in ([^ two doubly degenerate 
eigenvalues are found, of equal magnitude but opposite 
sign: ± y^i732CT^6+^ . The positive one is associated 
to ti = cos and the negative one to ti = sin. As for the 
cases (ni = 0, 712) and (ni, 712 = 0), in which /3+^ = and 
= respectively, the eigenvalues read a/tt/So^&q 6"^ 
and — ^/tt/Sct^&q 6"^, respectively. Finally, the eigenvalue 
associated to the fundamental Gaussian eigenvector, that 
is, to the case (ni = 0,n2 = 0), reads ^/^^/2a^bQbQ . 

We have therefore shown that one is able to master 
the number of supermodes and the magnitude and sign 
of their eigenvalues by a proper choice of the modulation 
amplitudes of the parametric multimode interaction ker- 
nel (PT|) . In Figure [T] we give two examples of kernels, 
designed so as to lead to eigenvalues respectively pro- 
portional to {1, 1, —1, —1} and {—4, 1, 1, 1, 1}. We note 
that the results shown in the figures have been obtained 
by numerically diagonalizing the kernels, what gives ad- 
ditional support to our previous approximate analytical 
treatment. 

Our analysis has been restricted so far to the case of 
kernels with symmetric Gaussians. However, it can be 
readily extended to the more physical case of a kernel 
factorizing in two Gaussian functions of variables x + x' 
and X ~ x' having unequal widths [13, . In fact this 
is the most general situation as it is not always possible 
to configure the OPO/OFA in such a symmetric way (in 
SPOPOs, for instance, for typical situations one has a 
difference of one order of magnitude between Gaussian 
widths). In this case the kernel we propose has again the 
factorizcd form now with 

K± (x) = c-^-±-' COS {P^x) , (24) 

instead of ((2T|) . where we allowed for different widths, 
(T^^, along the directions x + x' and x — x' . In Appendix 
B we give the eigenvalues and eigenvectors of this general 
case. 
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FIG. 1: (a) Density plot of the kernel — see (|19[) and (|2ip — particularized to the following parameters: a = 0.005, — 1, 
= Sttct, = 1, 6^ = 1 /3j = Gttct, the rest of b's being zero; in (b) we show its associated eigenvectors with the 
corresponding eigenvalues, (c) Density plot of the kernel for a different choice of the parameters: a — 0.005, = —2, b^ — 1, 



1, P-^ = Sttct, &2 = 1; = Bttct, the rest of 6's being zero; in (d) we show the corresponding eigenvectors and eigenvalues. 



IV. CONTROLLING THE NUMBER OF 
SUPERMODES AND THE EIGENVALUES IN 
ACTUAL PARAMETRIC DEVICES 



The method we sketched out above can be practically 
applied to spatial or temporal modes of the OPO/OPA 
since in both cases a factorized form (fT9)) of the kernel is 
a very good approximation (it is even exact in the "spa- 
tial case" ) , and both K± admit as well an approximated 
Gaussian form, see e.g. [13, IH, HI]. The point is then 



how to implement in real devices the quasiperiodic mod- 
ulations wc introduced in the previous Section. 

We will divide our presentation in two Subsections. 
The first one deals with what we call the continuous case, 
which includes both problems in which a continuous bo- 
son representation is used (single-pass devices) or prob- 
lems in which, being that representation discrete, the ker- 
nel can be treated as a continuous function (SPOPO [i^ 
and the self-imaging OPO The second Subsection 

deals with what we call the discrete case, in which a dis- 
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Crete boson representation is used and a discrete treat- 
ment of the problem is simpler, given the relatively small 
number of modes involved, like an OPO tuned to a single 
transverse mode family [37 1. 



A. The continuous case 

As for the Gaussian form of the subkernels K± wc note 
that both the phase-matching function D (|4I5|) of the 
temporal case and the diffraction function A ([7]) of the 
spatial case can be well approximated by the Gaussian 
e-5^?("+'^')'-5'^- In the temporal case [H,!!! 



Tl 



(25) 



/Cp g being the derivatives of the pump/signal wavenum- 
ber with respect to frequency at phase matching, while 
in the spatial case ti = 0. Hence if the pump spec- 
tral amplitude ap (w) has a Gaussian shape e^^'^p'^ ^ Tp 
being the individual pulse duration, the total kernel can 
be approximated hy K = e-5'^+('^+"')'-5'^-('^-'^')' with 
CT^ = Tl +Tp. If the pump has not a Gaussian shape but 
can be expressed as a sum over cos functions, Tp = 
(hence cr^ = ) and the Gaussian kernel will be multi- 
plied by those modulations (see below). 

As the pump amplitude ap is a factor of see Eq. 
([3]), a simple and practical way to tailor the kernel is by 
shaping the pump, as we will see more precisely now. 



diffraction gratings perform the inverse Fourier trans- 
form so that the outgoing field has the following tem- 
poral shape: ap,out(0 = M{t) (gD ap^in(i), i.e. the usual 
convolution, which is used to pump the SPOPO. It is 
then straightforward to see that, by using ap^out(w) in 
Eq. ([3]), one obtains a modulation of the kernel K{x,x') 
as in Eq. (|24p . with 6+ — &„, 6q = 1 and &~>o =0. As a 
consequence, the eigenvalues A„ are proportional to ±&„ 
according to the analytical solutions given in Section Iml 



b. Using pulse shapers leading to special waveforms: 
In this case we assume that a pulse shaper has been pro- 
grammed so as to convert the frequency comb coming 
from the pumping laser (whatever its temporal wave- 
form could be) into into a sequence of equal rectangu- 
lar pulses at the frequency comb repetition rate. Hence 
the pump spectrum is proportional to sine (Tpw), Tp be- 
ing the pulse duration. Then, if is much smaller 
than the width of the crystal response function along the 
direction iv + lo' — equal to Tf ^ defined in Eq. ([25|) — , 
i.e. for "long" pulses, the kernel can be approximated by 

e-i-i(-W)=-i^^(<^-<^')%inc[Tp(w + w')], see Eq. ©, 
where fT+ ~ ti [43 |. As a consequence, as sine (Tpw) can 
be approximated by a finite (Fourier) sum of functions 
cos(n^w), with L large as compared to t^^ , the ker- 
nel reads as in ([24| . As the coefficients of the Fourier 
series, 6^, are all very similar (because of the function 
sine (tw), whose Fourier transform is rectangular) we ex- 
pect a large degree of degeneracy between eigenmodes, 
as is actually evidenced by numerically diagonalizing the 
corresponding kernel. 



1, Tailoring the pump 

We focus on the temporal case psi . |43 . |49| and hence 
deal with pumps consisting of pulses. However the ideas 
put forward below can be applied equally to the spatial 
case by substituting the pulse shapers we consider 
by amplitude masks or simply by superposing different 
plane waves. 

a. Using pulse shapers with harmonic spectral re- 
sponse: A first possible method consists of using a pulse 
shaper [57| with several well separated amplitude modu- 
lations on the Spatial Light Modulator (SLM). The effect 
of such pulse shaper in the frequency domain can be rep- 
resented in the general case by the modulating function 

N 

A/(w) = ^&„cos(/3„Lj), (26) 

where /3o = by definition. The incoming pump field, 
of temporal envelope ap.in(i)i is sent onto a diffraction 
grating and the image ap.in(w) is created in the Fourier 
plane of a lens in correspondence of the SLM mask. Im- 
mediately after the mask, the outgoing field has a spec- 
trum ap,out(w) = M(a;)ap^in(w). Then another lens and 



In Figure [5^ we trace the kernel corresponding to a 
realistic situation where a SOPO cavity, with the nonlin- 
earity of a 100/^m-thick BIBO crystal (ti 20fs; see 
for details), is pumped by a train of square pulses. The 
eigenvalues A„ obtained by the numerical diagonalization 
of K for increasing values of Tp are shown in Figure [2]d 
and are compared with the corresponding cases of Gaus- 
sian pump pulses: as expected, the degree of degeneracy 
increases dramatically with the duration of the square 
pump pulses. For instance, for Ips pulses the first hun- 
dred eigenvalues differ by less than 1%, see the inset in 
Figure Oa, what allows the generation of very high qual- 
ity, highly multi-dimensional entanglement. 

c. Using delay lines A third method consists in us- 
ing delay lines that, starting from one comb, allow the 
superposition of several combs so that each "tooth" in 
the temporal domain is made of a series of pulses de- 
layed/advanced by tn with respect to the first one; one 
then gets as many eigenvalues as there are different super- 
posed combs, the eigenvalues being proportional to the 
amplitude of each comb as we show next. For example, 
a{t) being the (normalized) envelope of the pump field 
associated to a train of pulses incoming the OPA or the 
SPOPO cavity, the envelope of the pump field associated 
to a generic superposition of delayed / advanced trains of 
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FIG. 2: (a) kernel evaluated in the realistic case of a 100/xm thick BIBO crystal inside a 4m ring OPO cavity pumped by a 
c.w. train of 500fs rectangular pump pulses for a degenerate type I critically phase matching operation at 0.4/^m pumping, q 
and m denote the indices of the cavity longitudinal modes, (b) Eigenvalues obtained by numerical diagonalization of the kernel 
shown in (a) when different trains of pump pulses with Gaussian (grey-dashed lines) and rectangular (black-solid lines) shapes 
are considered; the duration of the pulses are specified in the figure. Ai refers to the eigenvalue with the largest absolute value, 
note that the caption is a zoom of the upper-left corner of the figure. 



pulses reads 

ap(t) = 6o a{t) +^hn[a{t~ i„) +a{t + <„)] , (27) 



N 



n=l 



where bo and fo„ are real coefRcients controllable in the 
experiment. In the Fourier domain Eq. (|27p reads 



ap{uj) = a{uj) 



N 



6o + ^ &„ cos(t„a;) 



n=l 



(28) 



a{uj) being the Fourier transform of a(t). Since we have 
the same eigenproblem as in the previous cases, the eigen- 
values A„ will be again proportional to ±6„. 



2. Tailoring the crystal response 

Another interesting possibility in the temporal case 
is to modify not the pump temporal shape but the ef- 
fect of the nonlinear medium, which allows us to con- 
trol the function It can be accomplished by using 
several identical non-linear crystals which are not per- 
fectly phase-matched and separated by fixed distances. 
In the general case where N crystals of thickness Ic are 
used, centered at planes z = Zn, the corresponding phase- 
matching function D, see reads 

D{io,io') ^ I dze'^H-'")^ (29) 



N 

^exp 



i^$(^,^') 



sine [$ {uj,uj')] , 



where 3> (w, w') is given in ([S]). If the crystals are arranged 
symmetrically by couples at distances d„ between mid- 
planes, then the above expression becomes 



N/2 

D{lo,lu') = 2^cos 



sinc[$(a;,a;')] , (30) 



where now the sum extends over couples of crystals"^. 
When the pump spectrum is narrow as compared with 
the width of D along the direction uj + uj' — the quantity 
defined in Eq. (j25p — a safe approximation consists 

in setting $(w,aj') — > ti{uj + uj') in cos y!-<I> (w, w') 

Hence, 



cos 



dn 



$ (w, LU') 



COS 



dn 



Tl {uj + Cj') 



acts as a harmonic modulation (equivalent to 6+ = 2, 
f^n — T^'^'ij ^0 ~ 1' ^'^"^ ^n>o = in our previous no- 
tation) of the single-crystal kernel, and the latter can be 
approximated by a factorized Gaussian form as already 
discussed, so the full kernel takes the form ((M)) . 



^ If the number of crystals is odd, say N = 2M + 1, then Eq. 
I I30I I is modified by substituting the sum upper limit by M, and 
adding a term equal to sine [<1> {ui, ui')] to the result, as can be 
easily checked. 
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Discrete case 



As commented in Section IIIBl the discrete represen- 
tation is useful when deahng a low number of signal 
modes. A typical configuration corresponding to this case 
is given by OPO with monochromatic pump, whose cav- 
ity is tuned at the subharmonic to some transverse mode 
family represented by the family index / [s^]- Recall 
that family / contains the / + 1 Laguerre-Gauss modes 
i(/-i)/2,±i(r) with I = Iq, lo + 2, /, being Iq equal to 
for even families and 1 for odd families. 

In order to simplify the upcoming discussion, we as- 
sume that the the pump beam has a coaxial cylindri- 
cal symmetry, in which case Orbital Angular Momentum 
(0AM) conservation ensures that the down-converted 
photons must have opposite OAMs. This implies that, if 
one uses the basis of Laguerre-Gauss modes TEMpj, the 
parametric down-conversion Hamiltonian takes the form 



XI 



H.c. 



(31) 



with 



Xi 27r 



+ 00 



rdrap(r) 



7^: 



(/-0/2 



(r) 



(32) 



ap(r) and ''^(/_;)/2(^) being the transverse profiles of 
the (normalized) pump field and the Laguerre-Gauss 
modes * at the cavity waist, respectively, and a[ the cre- 
ation operator associated to the Laguerre-Gauss mode 
i(/_|i|)/2,i(i")j which we abbreviate as Li(t) from now on 
because a fixed value of the family index / is assumed as 
explained. 

The continuous boson operators a(r) defined on the 
transverse plane can be expressed in terms of Laguerre- 
Gauss modal operators as 



a(r) 



E — 



Sis 



[Li{v)ai + L_i(r)a_i] 



(34) 



The kernel K(y, r') of Hamiltonian ([8]) has therefore the 
following expression: 



(35) 



Explicitly, we have 



2p\ 1 / v^i 



tt{p + 1)1 w \ w 



(33) 



where Lp{x) are the modified Laguerre polynomials and to is the 
spot size of the beam at the cavity waist. 



This expression is close to the one assumed in (|T9)) . It can 
actually be brought to that exact form by introducing the 
Hybrid Laguerre-Gauss modes 



Siir) 



V2(l+'5i,o) 
1 

1^2(1+5,, o) 



[Li{r)+L^i{r)], 



(36a) 



[Li{r)-L^i{r)], (36b) 



which are respectively proportional to cosZ0 and sinZi/). 
The kernel can then be written as 



i^(r,r')= J2 E^'^'W^'(^')- 



(37) 



H=C.S I 



In other words, we have just shown that the Hybrid 
Laguerre-Gauss modes Ci and Si are the supermodes for 
the present OPO configuration. In addition, both the 
"sine" and "cosine" type modes have the same eigen- 
value Khj ~ xi {H = C,S). Note that nothing pre- 
vents xi from being complex, say xi = IXil 6xp(i^/';), 
what means that the supermodes with different I can 
be squeezed along different directions of phase space, di- 
rection (tt -I- ) /2 in particular for modes {Hi{y)}h=c,s- 
As for the number of available supermodes N , it is given 
by the number of modes contained in family /, that is, 
N = f + l. 

Let us now consider the situation in which the pump 
has a Gaussian shape, that is, ap(r) ~ WgGp{r)^ where 



1 



cxp 



(38) 



is a TEMqo mode with spot size Wp = pWs, Wg being the 
spot size of the transverse modes at the signal frequency 
and p a scaling factor (the factor Ws in ap is included for 
dimensional reasons, as we defined ap as a dimensionless, 
normalized pump amplitude). It is obvious from ([32]) 
that xio > XU)+2 > ■■■ > Xf: in other words, the lower 
the 0AM, the better the signal modes overlap with the 
pump profile. However, for a large enough p, the pump 
profile is somehow seen as plane by all the signal modes, 
and the couplings xi become independent of I: 



Xl ^ ^TTA / 

P»l V TT W 



+ 00 



rdr 



p JO 



7^' 



(/-0/2 



(r) 



(39) 

Hence, by pumping with a wide Gaussian, one can bring 
all the couplings xi to a common value, at the expense 
of decreasing them, hence increasing the OPO threshold. 
Note however that, as discussed next, the desired effect is 
obtained even for moderate values of p, so that the OPO 
threshold stays within reasonable limits. 

In Fig. [3] we show the ratios r; = xi/ XUs ^ function 
of p for the first families. Note that the p axis starts at 
1 / v^, which is the value one has in the doubly resonant 
configuration (cavity resonant both for the signal and 
the pump). Note also how r; goes very quickly above 
0.5, which is the value above which one finds more than 
90% squeezing for the corresponding supermodes (37j . 
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In the insets, we show i?th(p) = xf^ip = / V2) / xf^ip) , 
which gives the ratio between the pump power needed to 
make the signal field oscillate for a given p and that for 
p = l/-\/2; note that the threshold is not dramatically 
increased for reasonable values of p (it is even decreased 
for small p). 

One can even tune the actual values of the couplings 
Xi to specific values by pumping not with a single TEMqo 
mode, but with a combination of TEMqo modes of differ- 
ent widths. To show this, let us consider two examples 
within the third family of transverse modes (/ = 3, and 
then I ~ 1, 3). First, we want to make xi = ~X3j so that 
the eigenvalues of the / + 1 = 4 supermodes have the 
same magnitude, but opposite sign. In this case, we can 
use the following pump profile 

ap{r) ^ Ws[Ga{r) cosO - Gbir) sine], (40) 

that is, a superposition of two Gaussians with opposite 
phase, and spot sizes Wa = PaWs and wt = ptWs- Note 
that this type of pump shape can be easily obtained in the 
lab by mixing the Gaussian beams on a beam splitter. It 
is then straightforward to show that choosing the mixing 
angle such that 

..„,.f£.)Vi±Myi±i4, ,41) 



one gets xi = ~X3- Note that, in the general case where 
the resonator is tuned to family /, one can have complete 
control of all the ratios between the coupling parameters 
by using 1 + (/ — lo)/2 TEMqo beams with adjustable 
amplitudes as pumping beams. 

In the second example wc show the very interesting 
case where some of the couplings become zero by using 
an appropriate pump shape. To show that this extreme 
case is indeed possible, we consider the previous example 
with the following choice of the mixing angle between the 
Gaussians: 

in which case it is simple to show that xi = Oi while X3 
can be large enough with a proper election of pa,f)- It can 
be interesting, for example, to choose the mixing angle 
between the Gaussians in order to cancel the coupling 
with the TEMqo modes in a general OPO, so as to favor 
the coupling to the TEMiq/qj modes, thus forcing it to 
emit the signal and idler modes with opposite 0AM. This 
way, one can induce a spontaneous breaking of the radial 
symmetry, what has been predicted to give rise to some 
remarkable quantum properties [35l |4]| . 

The methods discussed so far allow us to tune at will 
the eigenvalues of the supermodes, but not the spatial 
profile of the supermodes themselves (they are the Hy- 
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brid Laguerre-Gauss modes in all the cases). In order 
to change the form of the supermodes, we can go a lit- 
tle further and add beams with non-zero 0AM to the 
pump, what would allow us to engineer any kind of cou- 
pling between the modes with definite 0AM of a given 
family. For example, if one uses a pump field with a +2 
0AM component, it becomes possible to couple the sig- 
nal modes with +1 and —1 + 2 0AM; mixing several of 
these pump beams, one can even tune each of these cou- 
plings to a desired complex value, so much as we have 
shown for the couplings between the signal modes with 
±/ 0AM by pumping with zero 0AM Gaussian beams. 




FIG. 4: Examples of cluster states. The dashed line denotes 
that the corresponding coupling has opposite sign. 



V. APPLICATION TO THE GENERATION OF 
ARBITRARY CLUSTER STATES 

Having a set {sn{£.)}n=i,2,...,N of copropagating modes 
with squeezing properties chosen at will (as we have dis- 
cussed along the previous sections) offers the possibil- 
ity of generating any type of continuous- variable cluster 
state, that is, of producing any type of Gaussian multi- 
partite entangled state. To see this, just note that we 
have proved that these modes, which we have called su- 
permodes, evolve according to the Hamiltonian 



= E ^-Sl' + H.c. 



N 



(43) 



where is the creation operator associated to the su- 
permode s„{£,) of the kernel K{^,£^') having eigenvalue 
A„, both of which can be controlled via any of the ideas 
explained in the previous section. Note that this Hamil- 
tonian can be written as 



(44) 



where St (S-J , S*], S"]^) is a row vector operator, su- 
perscript "T" denotes transposition, and we have defined 
the diagonal matrix C = diag(Ai, A2, A^r). Through 
an arbitrary unitary matrix we can define a set of cre- 
ation operators = {bI, bI, b]^) ~ S^'U^ for some 
new modes, so that 



^i=i^Bt^(Bt)T+H.c., 



(45) 



where K. = UCU^ is a new coupling matrix. Hence (|43|) 
can be seen as the diagonal form of (|45p . For each choice 
of the coupling matrix /C, which can be chosen as sym- 
metric without loss of generality, this Hamiltonian gen- 
erates a different type of multipartite entanglement be- 
tween the B modes. 

As simple examples, consider the cluster states repre- 
sented in FigureSl In the first type, four modes of a circle 
are connected (all with same strength) only to their first- 
neighbors, one coupling having opposite sign respect to 
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the rest. The corresponding coupling matrix is (appart 
from a multiplicative factor) 



(46) 



which has doubly degenerate eigenvalues ±1. Note that 
these eigenvalues correspond exactly to some of the ex- 
amples considered in previous sections, in particular, to 
the example shown in Figure [T^ of Section IIIIl and to 
the second example of the previous section, where the 
couplings xi s-nd X3 of an OPO tuned tuned to the third 
family of transverse modes were tuned to the same abso- 
lute value and opposite sign. 

In the second example of Figure^ five modes are inter- 
connected with the same strength via the coupling matrix 
(again appart from a possible multiplicative factor) 



K. = - 
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(47) 



which has eigenvalues {—1, 1/4, 1/4, 1/4, 1/4} coinciding 
with the example of Figure [Ud in Section IIIIl In the dis- 
crete spatial case, these eigenvalues can be obtained by 
tuning the resonator to the / = 4 family of transverse 
modes and tuning the couplings to xo = —4x2 = —4x4, 
which can be done via the multi-Gaussian pump tech- 
nique discussed in the previous section, in this case with 
1 -I- (/ — lo)/2 = 3 Gaussians of appropriate widths and 
weights. 

Incidentally, we note that these two cluster states are 
of very different nature; while in the first case tracing 
out one of the modes does not destroy the bipartite en- 
tanglement shared between the rest of the modes, this is 
not the case in the second example, in which tracing out 
one of the modes leads to a completely separable state. 
In the latter case, we then say that the entanglement is 
"genuinely multipartite" 0. 

Indeed, the state of the second example in FigureJH 
corresponds to the class of so-called GHZ-like states [3], 
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FIG. 5: Transverse profile of tlie modes sliaring genuine mul- 
tipartite entanglement in the first families; they are obtained 
by applying a Braunstein rotation onto the Hybrid Laguerre- 
Gauss modes (which act as the supermodes in our system, 
that is, as independently squeezed modes). 

that is, states which tend to 

|GHZ)jv ^ J dx\ x,x,...,x) , (48) 

N 

in the limit of perfect squeezing (kind of generaliza- 
tion of EPR states to N modes), and are then com- 
pletely characterized by the variance of the (N + 1) joint 

quadratures Y.f=i ^3 |p, - Pj+i(modAr) | _^ , being 

Xj = b] + 4 and Pj = i [b] ^ Bj^ HI, of which 

|GHZ)7v is indeed an eigcnstatc. 

There is another way of obtaining GHZ-like states from 
a set of N squeezed modes I Sn I f^— i ; which was intro- 
duced by van Loock and Braunstein [lOj, who proposed 
to start from N — 1 modes squeezed along a given di- 
rection of phase space, plus another one squeezed along 
the orthogonal direction — a state that can be produced 
in the way explained in the present paper — , and to ap- 
ply to this state a specific transformation known as the 
so-called "Braunstein rotation" [l^, ■ 

Note that in order to have a reasonable approxima- 
tion of the GHZ state, one needs large squeezing in all 
the initial modes. In the previous section we actually 
showed that by pumping with a wide enough Gaussian 
mode an OPO tuned to the /-th family of transverse 
modes, one can bring together the thresholds of the / + 1 
modes contained in that family, and hence all the super- 
modes will leave the OPO highly squeezed when work- 
ing close to threshold. Note however that all the su- 
permodes arc squeezed in their P quadrature (taking the 
pump phase as the reference), and hence, before applying 



Braunstein's rotation, one needs to perform a 7r/2 phase 
shift in one of them to obtain the true entangled modes. 
In Figure Owe show the transverse profile (square mod- 
ulus) of these modes for the first families (the 7r/2 phase 
shift has been applied to the Cig supermode). Note that 
using SPOPOs pumped by rectangular pulses we showed 
that one can generate hundreds of highly squeezed su- 
permodes (sec Figure and hence, this system offers 
a highly dimensional alternative to the OPO tuned to 
a given transverse family where the number of squeezed 
supermodes is more modest. 



VI. CONCLUSIONS 

We have shown that an appropriate shaping of the 
pump beam (or of the geometry of the nonlinear medium) 
enables us to generate in a multimode OPO any Hamil- 
tonian bilinear in the annihilation operators of the down- 
converted modes, and therefore to generate highly multi- 
mode non-classical states of light that may be of interest. 
In addition, we have shown that these quantum states 
are easily characterized in a special mode basis, the ba- 
sis of "supermodes" , as a superposition of independent, 
squeezed copropagating modes. 

The measurement of the quadrature components of any 
mode is always possible with a balanced homodyne set- 
up using a coherent state in the mode of interest as a 
local oscillator, which projects the multimode state on 
the mode of the local oscillator. However, this technique 
is destructive, so that it allows for a measurement of the 
properties of one supermode at a time. An alternative 
to this single-mode homodyne detection is multiplexed 
homodyne detection [H, |60|, with which it is possible 
to recover simultaneous information of the different su- 
permodes using an appropriate data processing protocol. 
This opens the possibility of considering the copropa- 
gating entangled modes as a valid resource for one-way 
quantum computation, at least in the Gaussian domain. 

Even though having a cluster state between coprop- 
agating optical modes could still be useful for certain 
quantum information processing tasks, a highly relevant 
question which deserves further consideration is the way 
to access the entangled modes independently (for exam- 
ple to perform measurements on each of them, what is 
needed for one-way universal quantum computation), or 
to separate them to further use in quantum communica- 
tion networks. 

There are different possible techniques to separate the 
different modes. For example in the spatial domain, 
diffraction gratings are often used but they introduce 
losses, which are detrimental for the quantum effects; 
combinations of spatial light modulators are promising 
candidates on this regard, as they are essentially unitary 
transforms on the spatial modes [6l|. In the time do- 
main, the well-known pulse shaping techniques can be 
of interest, but cannot transfer energy between different 
frequency modes; on the other hand, light modulators 
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can do this task, and are the equivalent of diffraction 
gratings, but they cannot modify the hght spectrum by 
a large amount. Then, for temporal modes, the solu- 
tion proposed in [62| of using sum-frequency conversion 
with an appropriately shaped pump beam is certainly the 
most promising. 
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Appendix A 



expansion. In such case, by defining 

s (x) = /cZi;i {x) + /sZl'j (x) , (52) 

where /c and fs are constants to be determined, and com- 
puting F [s (x)] one obtains 

F[.s ix)] = b+b^[ik^j, + k:fs)z'i:i {x) 

cell +mc + klfs)zl-:A^)], (53) 
where fcg = k'^'^'i \ and so on, arc short notations for the 

constants introduced in (|51|) . Then the eigenvalue equa- 
tion F [s {x)] = As (x) is trivially fulfilled by demanding 

b+b- {klu + Kh) = A/c 1 
b+b^ {kif, + kiu) = A/, / ' 

which determines the eigenvalues A by imposing that the 
system of equations has nontrivial solutions. Two A's are 
obtained (Ai_2), proportional to the product and 
for each of them the ratio /s//c becomes fixed, which 
defines the true eigenvectors, si^2 hi this case. 



The set of functions defined by (|22p allows the diago- 
nalization of kernel under (pi]) . It is simple to show 
that the action of the kernel on such functions, defined 
by the Fredholm integral F p^ . is 

F (^)] = E E (^) , (49) 

m— n2— 

where c = cos, s = sin, 

C ^ \ _ I,C,C,Tli,Tl2 c,c / \ I jLS,s,ni,n2 ^s,s / \ 

■^mi,m2 V*^/ "'c,c,mi.m2 ni-K'jV'^y ""c.c.mi ,m2 rii ,n2 V"^/ ' 

(50) 

and 



7 ti,t2,ni,n2 
"'t3,t4,mi,m2 



xh{/3+,y)U{p-^^y), (51) 



are constants. When computing F [zf^_^ „^ (x)] 
a similar relation is obtained, showing that 

{^mi,m2 (2^) '^mi,m2 (2^)} ^O™!^ ^ ^loscd SCt frOm 

which actual eigenvectors can be found. Analogously 
the set {z^^ (x) ,zf^^ (x)} is closed with respect 
to the Fredholm integral P^ . 

The method can be easily visualized in the simplest 
case where bf are the only non-null coefficients of the 



Appendix B 

Although cumbersome it is straightforward to show 
that, if condition (|23l) holds, the eigenvectors of the ker- 
nel (HI]) with Eq. (IMl) are 



\i (/3+x) t2 {/3-^x) H,„(V2- 



(55) 

where ti=i,2 stands for any of the trigonometric func- 
tions cos or sin, m e N, r = y/a+a^, and H,„ (x) 
are Hermite polynomials. Their corresponding eigen- 
values read A^^^^*^^ ,„ = \rnb^_b~^ when ti = cos and 

A*l,t2 



^mbnj^n-^ whcu ti = siu, whcrc 



A777 - (-1) 



, ^^2 / a+ - a- 
2(<T+ +cr_) \a+ +a- 



(56) 



Hence also in this case one masters the number and eigen- 
values of the supermodes as well as their shapes. Note 
that for (T-|_ = (T_ = cr (symmetric Gaussians) Am/o = 
while Ao = ^7r/32cr2, and one recovers the results of the 
symmetric case. Note finally that when modulations are 
absent the previous eigenvectors and eigenvalues coincide 
with those for a SPOPO pumped by Gaussian pulses (4^ . 
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